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Abstract 

The scaled momentum distributions of charged particles in jets have been mea- 
sured for dijet photoproduction with the ZEUS detector at HERA using an 
integrated luminosity of 359 pb~^. The distributions are compared to predic- 
tions based on perturbative QCD carried out in the framework of the modified 
leading-logarithmic approximation (MLLA) and assuming local parton-hadron 
duality (LPHD). The universal MLLA scale, Aefj, and the LPHD parameter, k'^^, 
are extracted. 
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1 Introduction 



The formation of jets of hadrons can be described as a convolution of parton showering and 
hadronisation. Within perturbative QCD (pQCD), the parton shower can be described 
as long as the energy scale involved is sufficiently above the intrinsic scale of QCD, Aqcd- 
Hadronisation describes the process by which coloured partons become confined in colour- 
neutral hadrons. It cannot be described within pQCD. 

Perturbative QCD calculations can be performed using matrix elements up to a certain 
order in the strong coupling constant, a^. Alternatively, a resummation approach can 
be adopted, such as the modified leading-logarithmic approximation (MLLA) [1], where 
in addition to the fixed-order matrix elements, a subset of dominant terms of all orders 
in Us are included. In particular, pQCD based on the MLLA can be used to predict 
the multiplicity and momentum spectra of partons produced within cones centred on the 
initial parton direction. The MLLA may only be used to describe partons at scales above 
some minimum cutoff, Aeff > Aqcd- The value of Aes is predicted to be independent of 
the process considered. The local parton hadron duality (LPHD) [2] hypothesis predicts 
that charged-hadron distributions should be related to the predicted parton distributions 
by a constant normalisation scaling factor, n'^^. 

Tests of the MLLA have been performed before using data from e^e~ collisions at LEP [3, 
4] and PETRA [5], deep inelastic scattering (DIS) ep collisions at HERA [6,7], (anti-) 
neutrino-nucleon interactions from the NOMAD experiment [8] and pp collisions at the 
Tevatron [9]. In this analysis, the multiplicity and momentum spectra of charged hadrons 
within jets are studied using photoproduction (7^) in ep collisions, in which a quasi- 
real photon emitted from the incoming electron collides with a proton. The events were 
required to have two and only two reconstructed jets and the sample was enriched in 
events in which the photon interacted electromagnetically as a point-like particle. The 
analysis probes energy scales in the range 19 to 38 GeV, which spans the energy region 
between those accessed previously by the ZEUS, using ep DIS collisions [6,7], and CDF 
collaborations [9] . The quantities Aefr and n'^^ are extracted and their universality tested. 

2 The MLLA framework 

The MLLA includes all terms of order a" log^"'(-E'?\j and a" log^"~^(i?jn\t), where n is the 
set of positive integers and is the energy of the initial outgoing parton in the centre- 
of-mass frame of the incoming struck parton and exchanged photon. The "pi" superscript 
denotes a parton-level quantity. The MLLA accounts for colour-coherence effects between 
diagrams of the same order of by enforcing an angular-ordering scheme [10]. 
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The MLLA equations describe the momentum and multiphcity spectra of partons at a 
specified scale, Qq. They are only strictly valid for partons satisfying a;P' = <C 1, 

where p^^ is the 3-momentum of a parton in the centre-of-mass frame. For the MLLA 
predictions used here, the singularities were regularised by a single p^^'^' cut-off at scale 
Qq > Acff, where p^^^''^' is the transverse momentum with which the parton was emitted 
with respect to its parent. This is not the only possible way to regularise the MLLA; 
other forms lead to different predictions, particularly at low x^^ [11]. 

Predictions at the lowest valid scale, Qq = A^s, give the so-called limiting momentum spec- 
trum of partons, Z)'™'P' = where ^p' = ln(l/xP') and N^^ is the multiplicity of partons 
produced within a cone of opening angle, 6'p', measured with respect to the axis of the 
initial parton. The predictions assume that 6'p' is small. The shape of the predicted spec- 
trum is roughly Gaussian, although it falls rapidly to zero as —> In sin(6'P') / Acs 
a consequence of the regularisation scheme adopted. 

For a gluon-initiated parton-level jet, the function is written as [12] 



7r/2 

-Tv/2 



cosh a + (1 — 2() sinh a 



b sinh a 



B/2 



^?^F^^[cosha + (l-2C)sinha] 1 — , (1) 
h sinh a ^ ^ ^ M vr ' ^ 



where the parameters 6 = 9 and B = 101/81 are QCD constants defined in terms of the 
number of colours, ric = 3, and number of fiavours, Uf = 3. The symbols F and Is denote 
the Gamma and i?*^-order modified Bessel functions, respectively. The other variables 
are defined as 

F = ln(Efisin(^f)/Aeff), C = 1 - ^V^^, a = tanh-i(2C - 1) + ^r, (2) 

and, in the MLLA, -PhMlla = 1- The value of rij = 3 was chosen. The MLLA assumes 
massless partons and does not include any mass threshold effects for heavy fiavours. 
When a value of n/ larger than 3 is used instead, the theory is observed to give a poorer 
description of this and other data sets [12]. 

At leading order (LO), the peak position of the limiting momentum spectrum, ^p'^k? is 
predicted to be at 

?poak 

where c = 0.29. 



^p^eak = ^y + - C, (3) 



The limiting spectrum of quark jets, Dq™'^^ , is related to that of gluon jets according to 

r=)lim,pl _ ^ rSlim,pl / . \ 
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where r = -^gijct/^q-jet is the ratio of parton multiphcities in gluon- and quark- initiated 
jets. In the MLLA, r = Ca/Cp = 9/4 where Ca and Cp are the gluon and quark colour 
factors, respectively. 

Photoproduction samples contain both gluon- and quark-initiated jets, in the fractions 
denoted by eg and eq = 1 — eg, respectively. Thus, the limiting spectrum for partons in 
all jets can be parameterised as 

D^-'Pi = (eg + D^g. (5) 

Solutions to the MLLA evolution equations have also been made at so-called next-to- 
MLLA order. Each of these solutions partially accounts for orders not included in the 
equations above. With the next-to-MLLA corrections, -FhMlla and r differ from their 
MLLA values and both have a weak dependence on Ef^^^. Next-to-MLLA -Fhmlla and 
r values have been used in this analysis in the same way as they were by the CDF 
collaboration [9], wherein more details can be found. Their values were taken from three 
different next-to-MLLA calculations [13], which differ in the way the additional orders 
are accounted for, leading to some spread in the predicted -Phmlla and r values. Here, 
constant values of -FIimlla = 1-3 ± 0.2 and r = 1.6 ± 0.2 were used, with the theoretical 
uncertainties covering the spreads. 

The LPHD approximation relates the limiting momentum spectrum of partons to that of 
charged hadrons within jets, JJ^^'^''^^^ via 

^lim,ch _ ^ch^Um,pl _ ^ch _^ ^lim.pl _ K D'^^^^ (6) 

i.e. K = k'^^ (eg + (1 — eg)/r). Due to isospin invariance, k'^^ is expected to be approxi- 
mately 2/3. 



3 The analysis strategy 

To compare the parton-level MLLA predictions to measured hadron-level data, while 
assuming LPHD, each variable within the MLLA had to be estimated using a related 
hadron-level quantity. The hadron-level estimator for was chosen to be Ej^t = ^2j/2, 
where E-^ct is the energy of either hadron-level jet in the dijet centre-of-mass frame and 
M2j is the invariant dijet mass. The quantity p^' was estimated using the momenta of 
the charged hadrons, ptrk- The loss of the neutral hadrons is accounted for via the LPHD 
factor k"^. The MLLA variable Of was estimated using the opening angle of a cone 
measured with respect to the reconstructed jet axis, 6^. Accordingly, the quantity 
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given in Eq. 6, was estimated using the hadron-level multiplicity distribution of charged 
hadrons per jet, N^^, measured in bins of E'jet and in cones of varying 6c, differentially in 
^ = In (E'jet/IPtrkD- These dN^^/d^ distributions will be referred to as the ^ distributions. 

4 Experimental setup 

The data analysed here were collected using the ZEUS detector during the 2005 to 2007 
running periods, in which electrons^ were collided with protons with energies of Eg = 
27.5 GeV and Ep = 920 GeV, respectively, corresponding to a centre-of-mass energy, ^/s = 
318 GeV. The total sample corresponds to an integrated luminosity of 359 ± 9 pb~^. A 
detailed description of the ZEUS detector can be found elsewhere [14,15]. A brief outline 
of the components most relevant to this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [16], the microver- 
tex detector (MVD) [17] and the straw-tube tracker (STT) [18]. The CTD and MVD 
were operated in a magnetic field of 1.43 T provided by a thin superconducting solenoid. 
The CTD drift chamber covered the polar-angle^ region 15° < 6 < 164°. The MVD sili- 
con tracker consisted of a barrel (BMVD) and a forward (FMVD) section. The BMVD 
provided polar-angle coverage for tracks with three measurements from 30° to 150°. The 
FMVD extended the polar-angle coverage in the forward region to 7°. The STT covered 
the polar- angle region 5° < ^ < 25°. 

The high-resolution uranium-scintillator calorimeter (CAL) [19] consisted of three parts: 
the forward, the barrel and the rear calorimeters. Each part was subdivided transversely 
into towers and longitudinally into one electromagnetic and either one (in the rear) or two 
(in the barrel and forward) hadronic sections. The smallest subdivision of the calorimeter 
was called a cell. The CAL relative energy resolutions, as measured under test-beam 
conditions, were 0.18/a/E' for electrons and 0.35/ \/E for hadrons, with E in GeV. 

5 Event reconstruction 

A three-level trigger system was used to select events online [15,20,21]. At the first two 
levels, general characteristics of photoproduction collisions were required and background 



^ The word "electron" is used as a generic term for electrons and positrons. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. 
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from beam-gas events was rejected. At the third level, jets were reconstructed by applying 
the /ct cluster algorithm [22] to the CAL cells and a loose dijet selection was applied. 

In the offline analysis, the hadronic final state was reconstructed using energy-fiow ob- 
jects [23,24] (EFOs), which were formed from a combination of track and calorimeter 
information. This approach optimised the energy resolution and improved the one-to-one 
correspondence between the detector-level objects and the hadrons. The EFOs were cor- 
rected to account for energy losses in the dead material and were forced to be massless 
by setting the energy component equal to the magnitude of the three-momentum. 

Jets were reconstructed from EFOs using the kx cluster algorithm [22] in the longitudinally 
invariant inclusive mode [25] using the pt recombination scheme and with the R parameter 
set to i? = 1. 

Photoproduction events are characterised by the low virtuality, Q^, of the exchanged 
photon. At LO, photoproduction can be categorised as being either direct, if the photon 
interacts as a point-like particle, or resolved, if it fiuctuates into a partonic system, which 
then interacts with the proton. The LO direct photoproduction processes are boson 
gluon fusion, '-/g — > gg, and QCD Compton scattering, 7g — > qg. Important kinematic 
variables are the inelasticity, y, and the fraction of the photon momentum transferred to 
the hadronic final state, x^y. The variable be approximated using the observable 

x"^^, defined for a dijet event as 



ELi E?^'^ ex p(-r/j^^(^ 
" 2^ 



where E^r^^ and 77-''^* denote the jet transverse energy and pseudorapidity in the laboratory 
frame, respectively. A value of x^^ approaching one indicates an event from a direct-like 
photoproduction process. 



6 Event selection 

To remove non-photoproduction events it was required that: 

• the longitudinal position of the reconstructed vertex was in the range |^vtx| ^ 40 cm; 

• 0.2 < ?/jB < 0.85, where yjB is the Jacquet-Blondel estimator [26] of y; 

• no scattered electron was observed in the CAL with E'^ > 5GeVand ye < 0.85, where 
E'^ is the energy of the scattered electron and ye is the electron-method estimator of 
y [27]; 

• P^'"V/E^ < 2 GeV^/^ where P^'"" and Et are the reconstructed missing and total 
transverse momenta, respectively; 
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• I^CAL ~ ^calI < 6 ns, where I^cXl ~ ^calI is the difference between the arrival times of 
the first signals in the top and bottom halves of the CAL; 

• ^trk/^trk > 0.1, where -^t'rk/-^trk IS the ratio of the number of tracks fitted to the 
primary vertex to the total number of all tracks. 

To select an exclusive dijet sample enriched in direct events it was required that: 

• two jets were found such that: 

- the highest E^?* jet, labelled 1, had IV'^"] < 1 and > 17 GeV; 

- the second jet, labelled 2, had IV'^^^I < i and ^^.^^V^t" > 0-8; 

- the first and second jets satisfied — (j^'^^'^l > O.Qtt, where 0^°* denotes the 
azimuthal angle of the jet; 

• no third jet was found with |ryj'=*3| < 2.4 and E^r^^^ > 6 GeV; 

• xf' > 0.75. 

To ensure that the tracks were well reconstructed and not associated with secondary 
charged particles generated via nuclear interactions within the detector material it was 
required that: 

• the track transverse momentum was greater than 150 MeV; 

• the track pseudorapidity was between ±1.7; 

• the track passed through at least 3 CTD super layers; 

• the track was associated to the primary vertex. 

The requirement that there be two and only two jets roughly balancing in E^r^^ and in 
opposite hemispheres ensured that the events were LO-like, where the energy scale is 
well estimated using M2j/2. The x'^^ criterion was applied to minimise the influence of 
multi-parton interactions (MPIs) [28-30], which generate additional flnal-state hadrons 
and can disrupt the correspondence between the MLLA predictions and the data. After 
all the above selection, the data sample contained 23,449 events. 

7 Acceptance corrections 

Effects due to the limited detector and trigger acceptance, efficiency and resolution were 
corrected for in the data using a sample of events generated with the Pythia MC 
model [31]. The direct and resolved photopro duct ion processes were generated separately 
and combined in the ratio that best flt the x'^^ distribution in the data. The Pythia 
model includes the LO (2 2) matrix elements, approximates higher-order processes 
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using initial-state and final-state parton showers and simulates hadronisation using the 
Lund string model [32]. The CTEQ5L [33] and GRV-G LO [34] parameterisations were 
used to describe the proton and photon PDFs, respectively. The main sample included 
MPIs, simulated using the "simple model" [31] within Pythia, although the effects from 
MPIs were predicted to be negligible in the final sample. The detector simulation was 
based on Geant 3.21 [35] and included a complete simulation of the three-level trigger 
system. 

The data were corrected bin-by-bin to the hadron-level using factors extracted from the 
MC equal to the ratio of the predicted hadron- to detector-level cross sections. Here, 
the hadron-level was defined to contain all particles with an average lifetime greater than 
0.01 ns. The size of the bin-by-bin corrections were typically around 1.5. 

The normalisation of the ^ distributions was set such that the integral of the distribu- 
tions over the full ^ range equalled {N^^'), where {N^^") denotes the average hadron- 
level charged-particle multiplicity within jets, with the appropriate cone and energy scale 
criteria applied. The values of (-^e't) were extracted from the data by measuring the 
corresponding charged multiplicity distributions. These were corrected to the hadron- 
level using unfolding matrices derived from the Pythia MC sample. Full details of the 
procedure are described elsewhere [6,36]. 

8 Systematic uncertainties 

A detailed study [37] of the sources of systematic uncertainty associated with the measure- 
ment was performed. The dominant sources contributing to the systematic uncertainty 
on the ^ distributions are listed below (the numbers in parentheses refer to the maximum 
uncertainty observed in any one bin): 

• the ±3% uncertainty in the GAL energy scale, propagated to the ^ distributions by 
varying the GAL energies in the MG simulation accordingly (±4%); 

• the uncertainty simulating nuclear interactions in the detector material and the pro- 
duction of charged secondary particles. This was propagated to the ^ distributions 
by varying the difference between the number of tracks gained and lost due to such 
effects in the MG by a factor of 2 (±4%); 

• the uncertainty in the tracking efficiency, propagated to the distributions using the 
procedure described below (+5%). 

The MG slightly overestimated the number of tracks in the data, probably due to either 
the uncertainty in the hadronisation model or to inadequacies in the detector simulation. 
The unfolding procedure is only strongly sensitive to the detector-level simulation rather 
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than the hadron-level MC model and it was assumed that this was the sole cause of 
the excess. This systematic uncertainty was evaluated by randomly failing detector-level 
tracks in the MC with track rejection rates evaluated in bins of Ej^t, (^trk and 1 /ptrk? where 
6'trk is the polar angle between the track and the jet axis. The largest rejection rate was 
14%. The analysis was then repeated and the resulting difference in the ^ distributions 
was included in the systematic uncertainty. All the systematic uncertainties were added 
in quadrature. 

In the next section, several fits of the data are discussed. While nominally fitting the data 
and when evaluating the associated values, only the statistical uncertainties were con- 
sidered. The systematic uncertainties on the data were propagated, however, to the fitted 
parameters using the "offset method" . To apply the "offset method" , the fit is repeated 
for each source of systematic uncertainty, shifting the nominal data by the uncertainty 
attributed to that one source. The differences between the values of the parameters 
extracted from the nominal and the shifted data are then summed in quadrature and 
included as the total systematic uncertainty on the parameter itself. 

9 Results and discussion 

The ^ distributions were measured in five bins of E-^ct and in cones around the recon- 
structed jet axes with opening angles 6^ = {0.23,0.28,0.34}. The characteristic energy 
scales of the five Ej^t bins, E^^t = {19, 23, 28, 32, 38} GeV, were equated with the mean 
Ejct value for all events contributing to that bin. They are shown in Fig. 1. Each of the 
distributions are observed to be similar in shape and are roughly Gaussian with more 
pronounced upper tails. 

To assess the validity of the MLLA predictions using the measured ^ distributions, two 
approaches were adopted. The first, discussed in Section 9.1, was based solely on the 
position of the peak of the ^ distributions, ^peak- The second was based on the full shape 
of the ^ distributions and is discussed in Section 9.2. 

9.1 The ^peak analysis 

The values of ^peak were extracted from the ^ distributions using a three-parameter Gaus- 
sian fit. In accordance with previous analyses [6,9], the distributions were fit in the range 
fi^ ± 1, where fi^ is the arithmetic mean of the ^ distribution over the full ^ range. The 
explicit ranges and x^/dof values of the fits are given in Fig. 1. The x^/dof values range 
between 0.48 and 1.33 and hence indicate that the fits are reasonable. 
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Uncertainty in the ^peak values due to the choice of fitting range was added in quadrature 
to the total systematic uncertainty. It was evaluated by changing the fit range to fi^ ± 0.9 
and /i^ ± 1.1, leading maximally to a t^i'liYo systematic effect. The largest and only other 
source contributing more than 1% to the systematic uncertainty was the CAL energy 
scale, leading to a ^'^2 86^ effect. The extracted values of .^peak are given in Table 1 and 
are observed to increase as the energy scale or 6c increases. 

The ^poak values are shown in Fig. 2 as a function of fi sin {6^) , where the characteristic 
energy scale here is fi = E^ct- Also shown at their characteristic energy scales are data 
from the ZEUS ep DIS [6,9] analysis and the OPAL [3], TASSO [5], NOMAD [8] and 
CDF [9] collaborations. There is an approximately linear relationship between .^peak and 
In (.Ejet sin (6'c)). This relationship was tested by fitting the ^peak data, measured with 
9c = 0.23, with a straight line, parameterised as .^peak = ^ (ln(-E'jet sm{9c)))+B. In the case 
where only the new ZEUS 7p data were considered, the best fit values for the coefficients 
were found to be A = 0.56 ± 0.06(stat.)+°;°^(syst.)and B = 1.16 ± 0.09(stat.)l°:?^(syst.). 
The x^/dof of the fit was 0.51. 

A test of the same linear relationship was made using the global data set in Fig. 2. The best 
global fit values for the coefficients were found to be A = 0.682 ± 0.007(stat. © syst.) and 
B = 1.009 ± 0.019(stat. ©syst.), with a x^/dof of 0.77. Here, all systematic uncertainties 
were treated as uncorrelated. The globally-extracted parameters are consistent with those 
extracted from the ZEUS data alone. The ZEUS 7p points are systematically below the 
global-fit line, however the differences are within the total experimental uncertainty. 

The MLLA in fact predicts a small square-root correction to the perfect linear dependence, 
as seen in Eq. 3. Assuming Agg is constant within the range of energies probed, Eq. 3 
can be directly fit to the .^pcak data, treating Acs as a free parameter. In the case where 
only the ZEUS -yp data with 6c = 0.23 were considered, the best fit value was found 
to be Acff = 275 ± 4(stat.)l8(syst.) MeV. The xV^of of the fit was 0.70, indicating 
a good fit. When the global data set was considered, the best fit value was found to 
be Acff = 246 ± 3(stat. © syst.)MeV. In the global fit, all uncertainties were treated 
as uncorrelated. The x^/dof of the fit, with this simplistic error treatment, was 2.2, 
indicating some discrepancy. The globally extracted value of Aeg is not consistent with 
that extracted from the ZEUS data alone. 

The energy dependence of Aefr was studied by using Eq. 3 to map each ^peak value to a 
corresponding value of Aeff. The results, given in Table 2 and shown in Fig. 3 as a function 
of Ejct, show no evidence that A^s is dependent on the energy scale. A weak dependence 
was observed in the CDF data [9], which span a wider range of energy scales. However, 
the data do suggest that the value of Agfj is weakly dependent on 9c- Specifically, Fig. 3 
shows that the values of A^s extracted from the wider cone data tend to be systematically 
larger. This behaviour was also observed by the CDF collaboration [9]. Both the 6c and 
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Ejet dependence seen by CDF would contribute to the discrepancy observed when fitting 
Eq. 3 to the global data set. 

In Fig. 4, the values of Aefr extracted using the .^pcak data are shown as a function of the 
energy scale and compared to the previous results from ZEUS [6] using ep DIS collisions, 
and the OPAL [3], L3 [4] and CDF [9] collaborations. The values are all largely consis- 
tent in the energy scale region shown, supporting the prediction that A^s is a universal 
parameter. 

9.2 The ^-shape analysis 

The ^ distributions were also fitted using the predicted limiting spectrum, according to 
Eq. 6. The quantities K and Aeg were treated as free parameters during the fit. The fitted 
MLLA functions are shown in Fig. 5. The fits were restricted to the ranges indicated by 
the vertical lines and the x^/dof values of the fits are also given and lie between 0.34 
and 2.72. Typically, in each Ejet bin, the x^/dof increases as 9c does. The x^/dof values 
indicate that, while the theory does describe many of the features of the data in the 
fitting ranges, there are differences. Specifically, the rising edges of the ^ peaks are well 
described. However, the upper tails of the distributions are not adequately reproduced. 
The same was observed in e+e~ [3,4] and ep DIS [6] data and to a lesser extent in high-i?jet 
pp data [9]. This is likely due to the specific MLLA regularisation scheme used here and 
in the other aforementioned analyses. 

As discussed in Section 2, the MLLA regularisation scheme used here causes the partons 
to be cut-off at p^''^^ = Acfj, whereas the hadrons in the data are not. This leads to 
an intrinsic discrepancy between data and theory. The discrepancy is present for all 
^ > 0, however the magnitude of the effect is small at low ^ and increases until, for all 
^ > In (E'jet sin(^c)/AefT), there are only hadrons and no partons. 

A consequence of this discrepancy is that, in order to fit the data using Eq. 6, a relatively 
arbitrary upper fitting bound, had to be chosen for each ^ distribution. The criteria 
used to set ^+ were that the resulting fits were reasonably stable and that ^poak <^ ^+ < 
In (i?jet/250 MeV) was satisfied, where 250 MeV roughly corresponds to the values of Aeg 
extracted from the ^peak data. The finite experimental ^ binning was also a consideration. 
It was chosen to use ^+ = w^pcak + (1 ~ w) In (i?jet/250 MeV), with w = 0.25 for the 
nominal fits. The sensitivity of K and Acs to the choice of the fitting range was treated 
as a systematic uncertainty and was evaluated by varying w by ±0.1. This source of 
uncertainty strongly dominates the overall uncertainty on Acs, leading to a l}o^6% effect, 
although K was found to be largely insensitive to it. The same lower fitting bound, 
^_ = ln(2), was used is all cases and both K and Acs were observed to be insensitive to 
a variation of by ±15%. 
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The values of Agg extracted from the MLLA fits are given in Table 2. The results are 
in reasonable agreement with those extracted from the .^peak data, although the values 
extracted using the MLLA fit have larger uncertainties. The value of Aes from the MLLA 
method with 9^ = 0.23 and averaged over i?jct, weighting each data point based only on 
its statistical precision, is Acff = 304 ± 6(stat.)+^2(syst.) MeV. 

Values of Hch were extracted from the fitted K values using Eq. 6 and the values of eg 
predicted for each i^jct bin by the Pythia model. The eg values were roughly constant 
in ii^jct, at eg ~ 0.2. The Kch values are given in Table 3 and are shown in Fig. 6. The 
total uncertainty is dominated by the theoretical uncertainty associated with the next-to- 
MLLA correction factors. The data suggest a weak dependence on 9c. Specifically, as 
9c increases, so too does the central value of K^h- This is significant when the high degree 
of statistical correlation between the three 9c samples and the bin-to-bin correlation in 
the systematic and theoretical uncertainties are taken into consideration. The same is 
true for the Kch values reported by the CDF collaboration [9] , which were obtained using 
a different extraction method. The ZEUS data in Fig. 6 do not provide any evidence that 
Kch is dependent on Ej^t in the range probed. 

The value of Hch, measured with 9c = 0.23 and averaged over Ej^t, weighting the data 
points based on their statistical precision, was Kch = O.55±0.01(stat.)lg:|]^(syst.)l[5:^^(theo.) 
The Kch value extracted here is in good agreement with that reported by the CDF collab- 
oration, Kch = 0.56 ± 0.05(stat.)±0.09(syst.). To compare to the values extracted using 
e"'"e~ and ep DIS data and assuming no contamination from gluon jets, the values have 
to be scaled by rCp/EnULLACA ~ 0.55. This leads to values of Kch ~ 0.7. These other 
results were found with 9c effectively set to 7r/2 however. 

10 Summary 

The multiplicity distributions of charged particles within cones centred on jets have been 
measured as a function of = In {1/xp), where Xp is the fraction of the jet's momentum 
carried by the charged particle. These ^ distributions have been measured in five bins of 
i?jet and with three different cone opening angles, 9c, for 7p events containing two and 
only two jets, using 359 pb~^ of ep data. 

The peak positions of the ^ distributions, ^peak, were extracted and observed to increase 
roughly linearly with In (i?jet sin (^c))- A single value of intrinsic MLLA scale, Aeg, was 
extracted by fitting the .^peak data according to the predicted relationship between ^peak and 
In {Ejct sin {9c) /Acg). The best fit value was found to be Aes = 275±4(stat.)l8(syst.) MeV. 

The Ejct and 9c dependences of Acfr were studied by calculating a value of Agg from each 
■Cpcak data point. The value of A^s weakly depends on 9c but no E^^t dependence was 
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observed. The Aeff data are consistent with previously pubhshed data sets using different 
initial states, supporting the prediction that Acs is universal. 

The ^ distributions were also fitted using the limited momentum spectra predicted by the 
MLLA and assuming LPHD, in the regions where they are applicable. The theory largely 
described the data in these regions. The fitted MLLA functions were used to extract the 
value of Aeflf as a function of E^ct and 6c. The value extracted using this method with 
9c = 0.23 and averaged over Ej^t, was A^s = 304 ± 6(stat.)t32(syst.) MeV. 

The value of the LPHD parameter Kch was extracted as a function of Ejct and 6c from 
the fitted limited momentum spectra. Corrections based on next-to-MLLA theory were 
included. The value extracted with 6c = 0.23 and averaged over Ej^t, was = 0.55 ± 
0.01(stat.)t[!:[]^(syst.)l[5:J^(theo.). The value of Kch has a weak dependence on 6c and is 
consistent with the results published by the CDF collaboration. The data support the 
assumption that is universal. 
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Ejet ( GeV) 




^peak 


stat. 


syst. 




0.23 


1.99 


±0.01 


+0.02 
-0.02 


19 


0.28 


2.10 


±0.01 


+0.01 
-0.01 




0.34 


2.20 


±0.01 


+0.01 
-0.01 




0.23 


2.11 


±0.02 


+0.02 
-0.01 


23 


0.28 


2.21 


±0.02 


+0.02 
-0.01 




0.34 


2.32 


±0.02 


+0.02 
-0.01 




0.23 


2.22 


±0.04 


+0.03 
-0.02 


28 


0.28 


2.34 


±0.03 


+0.02 
-0.02 




0.34 


2.44 


±0.04 


+0.04 
-0.01 




0.23 


2.25 


±0.07 


+0.09 
-0.05 


32 


0.28 


2.36 


±0.06 


+0.10 
-0.03 




0.34 


2.56 


±0.06 


+0.07 
-0.05 




0.23 


2.40 


±0.05 


+0.04 
-0.08 


38 


0.28 


2.50 


±0.08 


+0.07 
-0.18 




0.34 


2.59 


±0.07 


+0.08 
-0.15 



Table 1: .^pcak values in the five E^^t bins using the three 9c values. The statistical 
and systematic uncertainties are also given. 
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^peak analysis 


^ shape analysis 


^jct (GeV) 


Be 


Aeff ( MeV) 


stat. 


syst. 


Aeff (MeV) 


stat. 


syst. 




0.23 


272 


±5 


+6 
-8 


304 


±4 


+7 
-32 


19 


0.28 


280 


±4 


+5 
-5 


298 


±4 


+21 
-25 




0.34 


289 


±4 


+6 
-5 


303 


±3 


+15 
-30 




0.23 


280 


±7 


+6 

-7 


307 


±6 


+10 

-32 


23 


0.28 


291 


±9 


+3 
-11 


305 


±6 


+23 
-32 




0.34 


297 


±8 


+3 

-9 


301 


±5 


+26 
-29 




0.23 


279 


±16 


+8 
-11 


285 


±12 


+8 
-19 


28 


0.28 


282 


±14 


+8 

-9 


294 


±10 


+7 
-29 




0.34 


292 


±17 


+5 
-17 


287 


±9 


+29 
-23 




0.23 


310 


±33 


+22 
-41 


298 


±15 


+25 
-40 


32 


0.28 


321 


±29 


+14 

-49 


302 


±13 


+26 
-41 




0.34 


283 


±24 


+21 
-28 


286 


±14 


+28 
-27 




0.23 


290 


±23 


+38 
-16 


311 


±15 


+13 

-52 


38 


0.28 


301 


±37 


+48 
-33 


287 


±21 


+42 
-32 




0.34 


319 


±36 


+31 
-38 


297 


±17 


+21 
-42 



Table 2: Aeg extracted at the five E-^^t points using the three 9c values obtained from 
both the ^peak and ^ shape analyses. The statistical and systematic uncertainties 
are also given. 
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Ejet ( GeV) 


0c 


.ch 


stat. 


syst. 


theo. 




0.23 


0.54 


±0.01 


+0.03 
-0.02 


+0.11 
-0.09 


19 


0.28 


0.59 


±0.01 


+0.03 
-0.01 


+0.12 
-0.10 




0.34 


0.63 


±0.01 


+0.03 
-0.02 


+0.12 
-0.10 




0.23 


0.56 


±0.01 


+0.03 
-0.02 


+0.11 
-0.09 


23 


0.28 


0.60 


±0.01 


+0.04 
-0.02 


+0.12 
-0.10 




0.34 


0.63 


±0.01 


+0.04 
-0.02 


+0.13 
-0.10 




0.23 


0.55 


±0.01 


+0.04 
-0.01 


+0.11 
-0.09 


28 


0.28 


0.59 


±0.01 


+0.04 
-0.04 


+0.11 
-0.09 




0.34 


0.61 


±0.01 


+0.04 
-0.02 


+0.12 
-0.10 




0.23 


0.56 


±0.02 


+0.04 
-0.04 


+0.11 
-0.09 


32 


0.28 


0.59 


±0.02 


+0.04 
-0.04 


+0.11 
-0.09 




0.34 


0.61 


±0.02 


+0.04 
-0.03 


+0.12 
-0.10 




0.23 


0.56 


±0.03 


+0.05 
-0.06 


+0.11 
-0.09 


38 


0.28 


0.58 


±0.03 


+0.04 
-0.04 


+0.11 
-0.09 




0.34 


0.61 


±0.03 


+0.03 
-0.05 


+0.12 
-0.10 



Table 3: k values extracted at the five E^ct points using the three 6c values. The 
statistical, systematic and theoretical uncertainties are also given. 
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ZEUS 




Figure 1: The ^ distributions in the five Ej^t bins using the three 9c values. 
The ZEUS data are shown by the solid squares. The inner error bars represent the 
statistical uncertainty. The outer error bars represent the statistical plus systematic 
uncertainties added in quadrature. Gaussian functions (solid line) have been fitted 
to the data within the regions indicated (dashed lines). The x^/dof of each fit is 
given on the plot. 
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Figure 2: ^peak a function of fisin{6c) , where fi denotes the characteristic 
energy scale for each specific process. The ZEUS '-fp data (solid circles) are shown 
along with ep data from the ZEUS collaboration (diamonds) and results reported by 
the OPAL (crosses), TASSO (triangles), NOMAD (stars) and CDF (open circles) 
collaborations. The inner error bars on the ZEUS points represent the statistical 
uncertainty. The outer error bars represent the statistical plus systematic uncer- 
tainties added in quadrature for all data sets. The data have been fitted with a 
straight line. 
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Figure 3: A^s extracted at the five Ej^t points using the three 9c values. The ZEUS 
data are shown by the solid points. The inner error bars represent the statistical 
uncertainty. The outer error bars represent the statistical plus systematic uncer- 
tainties added in quadrature. The points have been shifted horizontally for clarity. 
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Figure 4: A^s as a function of fi, where fi denotes the characteristic energy 
scale for each specific process. The ZEUS •jp data are shown by the solid circles. 
Also shown are ep data from the ZEUS collaboration and results reported by the 
OPAL, L3 and CDF collaborations. The inner error bars on the ZEUS 7p points 
represent the statistical uncertainty. The outer error bars represent the statistical 
plus systematic uncertainties added in quadrature for all data sets. 
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Figure 5: The C, distributions in the five bins using the three 6c values. 
The ZEUS data are shown by the solid squares. The inner error bars represent the 
statistical uncertainty. The outer error bars represent the statistical plus systematic 
uncertainties added in quadrature. The limited momentum spectrum predicted by 
the MLLA (solid line) has been fitted to the data within the regions indicated (dashed 
lines). The x^/dof of each fit is given on the plot. 
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Figure 6: k extracted at the five Ejct points using the three 6c values. The ZEUS 
data are shown by the solid points. The inner error bars represent the statistical 
uncertainty. The outer error bars represent the statistical, systematic and theoreti- 
cal uncertainties added in quadrature. The points have been shifted horizontally for 
clarity. 
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